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УДК 519.1 
Н.Н. ВОДОЛАЗОВ 


ОБ ОСОБЕННОСТЯХ ПОТОКА В СЕТЯХ 
С БАРЬЕРНОИ ДОСТИЖИМОСТЬЮ 


Рассматривается задача о нахождении максимального потока в сетях с барьерной 
достижимостью при помощи модифицированного алгоритма Эдмондса-Карпа. Осо- 
бенность предложенного алгоритма в том, что он позволяет находить поток в се- 
тях, для которых разработанный ранее алгоритм [1] не находил решения. 
Ключевые слова: сеть, источник, сток, поток, блокирующий поток. 


Введение. В классической постановке задачи стационарным потоком на 
графе С(№,И,Ф) называется отображение Х:И - К* , удовлетворя- 
ющее условиям: 


У; лад = У ла. Мое Г, + 5,0 #1, 


ие ®_ (5) ие ® , (©) 
(1) 
0< Л(ым) < с(м), 
где ® (0) - множество дуг, входящих в вершину ® ; 9 (0) - множе- 
ство дуг, выходящих из вершины ® ; с(м) - пропускная способ- 
ность дуги; 5 -— источник; 2 - сток. 
Требуется максимизировать > ла) ь 
иЕ®_ (5) 

Для решения этой задачи может быть применен симплекс-метод. 
Однако в 1956 г. Форд и Фалкерсон предложили более эффективный алго- 
ритм решения задачи, основанный на поиске увеличивающих цепей из ис- 
точника в сток. Эдмондс и Карп [2] в 1972 г. улучшили алгоритм Форда- 
Фалкерсона и получили алгоритм, время работы которого не зависело от 
пропускных способностей дуг. 

Постановка задачи. Требуется найти максимальный поток в графе 
С(И,И,Ф) с дополнительным ограничением на пути частиц, которое в 
[1] названо барьерной достижимостью. Множество дуг графа (/ разбито 
на три непересекающихся подмножества: С/,, (нейтральные дуги), С/, 
(увеличивающие дуги) и С/„ (барьерные дуги). С каждой частицей, про- 
ходящей по графу, связано целое число 7 - её барьерный уровень. Для 
частиц, выходящих из источника, он равен 0. Если барьерный уровень 
меньше некоторого заданного числа А , то при прохождении дуг из мно- 
жества (7, барьерный уровень увеличивается на 1, при прохо- 


ждении дуг из множеств С/„, и С/»ь барьерный уровень не изменяется. 


Частица может проходить по дугам множества С 5 только в случае, если 


её барьерный уровень равен А. 

Существующее решение. Решение описанной задачи предложено в [1]. 
Для этого по исходному графу строится вспомогательный граф. Каждой 
вершине Х исходного графа ставится в соответствие А+ 1 вершина 


127 


Раздел «Управление, вычислительная техника и информатика» 








вспомогательного графа х‘”, х“_...., х“?. Каждой нейтральной дуге 
(х, у) исходного графа ставится в соответствие А+ 1 дуга вспомога- 
тельного графа вида (х“”, у”), = 0....,А (рис.1). 


х, О } ЗУ» 
хьФ РОУ, 1 
мо > 
Хо ТО у, 


Рис.1. Преобразование нейтральных дуг 


Каждой увеличивающей дуге исходного графа (х, у») ставится в 
соответствие А+ 1 дуга вспомогательного графа: А дуг вида 


(х®, уб+0), {= 0,...,А-1 и одна дуга (х“, у“) (рис.2). 


Х, Ух 
Х к - р ЗУ: 


х ыы 
1 У! 
х 
о о Уо 
Рис.2. Преобразование увеличивающих дуг 


Каждой барьерной дуге исходного графа ставится в соответствие 
одна дуга вспомогательного графа (х“?, у“”) (рис.3). 


х, О м» 
х,_Ф — ЗУ, 
хо Ом, 


Хоо Ом. 


Рис.3. Преобразование барьерных дуг 


В соответствии с идеей, изложенной в [1], пропускная способность 
дуг вспомогательного графа не определена. Вместо этого определен макси- 
мальный поток по множеству дуг вспомогательного графа, который постав- 
лен в соответствие дуге исходного графа (множеству связанных дуг). Сле- 
довательно, применить алгоритм Форда-Фалкерсона мы не можем. 

В работе [1] авторами был развит новый подход к построению алго- 
ритма. Его сущность состоит в том, что алгоритм ищет кратчайший простой 
путь на вспомогательном графе из источника в сток, по которому можно 
пустить увеличивающий поток. Поток увеличивается на некоторую величи- 
ну, и пропускная способность множеств связанных дуг, через которые про- 
шел увеличивающий путь, уменьшается на эту же величину. Таким об- 
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разом, в отличие от алгоритма Форда-Фалкерсона алгоритм, описанный в 
[1], ищет не увеличивающую цепь, а увеличивающий путь. Ясно, что 
найденный поток будет неулучшаемым. 

Приведенный ниже пример показывает, что неулучшаемый поток 
может быть не максимальным в случае, когда возможен блокирующий по- 
ток [3]. Блокирующий поток в [3] рассматривается как поток, для которого 
не существует увеличивающего пути из источника в сток. Надо сказать, в 
таких случаях (блокировки) алгоритм Форда-Фалкерсона находит выход, 
поскольку существует увеличивающая цепь. Рассмотрим пример. 

Пример 1. Пусть дан граф С’ (рис.4). 





Я пе оп 


Рис.4. Нулевой поток на графе С 


На рис. 4 источник — вершина „5, сток — вершина 2. Первая цифра 
в записи О/1 - поток по дуге, вторая — пропускная способность. Все дуги 
нейтральные и нет барьерных дуг, т.е. мы рассматриваем случай поиска 
потока в обычном графе. Тогда первый увеличивающий путь, найденный 
алгоритмом, будет проходить по вершинам 5, <,./`.2 (рис.5). 





Я ше пл 


Рис.5. Блокирующий поток на графе С 


Поток на рис. 5 — блокирующий, так как нет увеличивающего пути 
из источника в сток. Величина потока равна 1, максимальный поток, изоб- 
раженный на рис. 6, равен 2. 
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ше ши! 


Рис.6. Максимальный поток на графе С 


Таким образом алгоритм не находит максимальный поток для реше- 
ния классической задачи. В случае с добавлением дополнительных ограни- 
чений на поток также возможно появление блокирующего потока. 
Предлагаемый алгоритм решения задачи. Для удобства будем 
рассматривать не суммарный поток частиц всех барьерных уровней через 
дугу, а поток частиц каждого барьерного уровня через дугу. Поэтому для 
поставленной задачи полагаем, что поток — — отображение 
Х:Ц > (В )^*', те. поток ставит каждой дуге в соответствие вектор 
чисел, при этом 2-я компонента вектора является потоком 2 -го барьер- 


ного уровня по дуге. Обозначим Л;(м,®) 2-ю компоненту вектора пото- 


А 
ка по дуге (м,0),а Л(и,0)= ). (и). 
1=0 


Условия (1) заменим на следующие условия: 


у. Ло (м) = ь» Лом), УюЕГ, + 5,0 #1, 


ие ®_ (пм ие ® , (0) 
Улад+ У Льард= У дао 
ие ®_ (пн ие ®_ (фи, ие ®_ (0) ь 


Г=1,... А-1, УОЕЙ; 


У лоад+ У лнад+ У Ллад+ У Лад= У, Ла 
ие ®_ (о )мИн ие ®_(0)пО, ие ®_(6)пИ ие ®_($)мИв ие ®_ (0) 
УшЕ И; (2) 
@а0>0, УЕ = 1.....А-1МиЕИ; 
О0< Л(м) < с(м). 
Требуется максимизировать уз Л (и)- №: У (4) ‚ т.е. вме- 
ие © _ (5) ие ©, (5) 
сто одного уравнения для каждой вершины в классическом случае, для 
данной задачи получается А - 1 уравнение для каждой вершины. Следует 
подчеркнуть, что эту задачу так же, как и классическую, можно решать 
симплекс-методом, но мы попытаемся изменить алгоритм Эдмондса-Карпа 
так, чтобы с его помощью можно было решать задачу о максимальном по- 
токе на графе с дополнительными ограничениями на пути частиц. Модифи- 
цировав алгоритм Форда-Фалкерсона, Эдмондс и Карп предложили произ- 
водить увеличение потока по кратчайшей увеличивающей цепи (считая, 
что все дуги имеют единичную длину) [2]. 
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Для решения поставленной задачи предполагается помечать каж- 
дую вершину вектором из А + 1 метки, где А - максимальный барьерный 
уровень; каждая метка содержит информацию об увеличивающей цепи со- 
ответствующего барьерного уровня в данную вершину. Наш алгоритм на 
первом этапе решения использует в своей работе очередь помеченных 
вершин и вспомогательный алгоритм пометки вершин. 

Алгоритм пометки вершин 

На входе помеченная вершина Х графа С’. Алгоритм помечает 
смежные вершины и добавляет их в очередь помеченных вершин. 

1. Для каждой дуги (х.а) и (а,х). 

2. Для каждой метки (В,с,“) 1-го уровня вершины Х. 

а. Если мы рассматриваем дугу Сх, “): 
1. если /(х,а)< СС, а) , то обозначим 
е = пищ(а, С(х,а)- Л(х,а)), 5='"+", иначе пере- 
хо- 
дим на шаг 2 и рассматриваем следующую метку. 
Ь. Если мы рассматриваем дугу (ах): 
1. если Л.(а,х)> 0, то обозначим 
е = пип(а, /,(х,а)), 5="-", иначе переходим на 
шаг 2 и рассматриваем следующую метку. 
с. Если (Хх, а) или (а, х) - нейтральная дуга и @ не поме- 
чена на 2 -ом уровне, то: 
1.  помечаем вершину @ на Ё-ом уровне меткой 
(х,5,е); 
П. если @ не первая в очереди, то добавляем @ в оче- 
редь помеченных вершин. 
Ч. Если (>, “&) - увеличивающая дуга, #< К и Я не поме- 
чена на (2+ 1) -ом уровне, то: 
1. помечаем вершину @ на + Ю-ом уровне меткой 
(5. а) 
П. если Я не первая в очереди, то добавляем @ в оче- 
редь помеченных вершин. 
е. Если (@,>х) - увеличивающая дуга, #> 0 и @ не поме- 
чена на (7 — 1) -ом уровне, то: 
1. помечаем вершину @ на @- Ю-ом уровне меткой 
(х,5,е); 
П. если @ не первая в очереди, то добавляем @ в оче- 
редь помеченных вершин. 
Ё. Если (х,а) или (а, х) - барьерная дуга, 1 = К и @ не 
помечена на А -ом уровне, то: 
1.  помечаем вершину Я на А-ом уровне меткой 
(х,5,е); 
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П. если Я не первая в очереди, то добавляем @ в оче- 
редь помеченных вершин. 
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Алгоритм поиска максимального потока 
1. Поиск увеличивающей цепи. 


а. Помечаем источник вектором | ® ® ® ‚ Прочерк на Е - 








ом месте означает, что вершина не помечена на данном ба- 
рьерном уровне. 


Помещаем источник в очередь помеченных вершин. 


Пока очередь просмотренных вершин не пуста: 


1. 
п. 


Ш. 


выбираем одну вершину Х из очереди; 

выполняем для нее алгоритм пометки вершин, если сток 
оказался помечен, то найдена увеличивающая цепь, 
переходим к шагу 2, иначе переходим к шагу 1; 

если очередь просмотренных вершин пуста, то алгоритм 
завершает свою работу. 


я, Поиск величины, на которую можно увеличить поток. 


а. 


Восстанавливаем увеличивающую цепь в сток по меткам. 
При этом цепь может проходить по некоторым дугам 
несколько раз, но на разных барьерных уровнях. 


Обозначим А величину © из метки стока. 


Проходим по полученной цепи. 


1. 


Пусть 7 - количество проходов цепи по дуге, при кото- 

ром направления цепи и дуги совпадают, Р - количе- 

ство проходов, при котором направления противопо- 

С(х,а)- (а) 
1-Р 





ложны. Тогда, если /> ри 4> 


С(х,а)- Г(х,а) 
1[-р й 


7! 





то присвоим # = 


3. Увеличение потока. Проходим по всем дугам цепи: 


а. Если направление цепи из источника в сток совпадает с 
направлением дуги, то увеличиваем поток по дуге на ве- 
личину д. 
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Ь. Если направление дуги не совпадает с направлением цепи 
из источника в сток, то уменьшаем поток по дуге на величи- 
ну д. 

Отличие предлагаемого нами алгоритма от алгоритма Форда-Фал- 
керсона в том, что, во-первых, он работает с вектором меток, а во-вторых, 
требуется шаг 2 для определения максимального потока по цепи. Необхо- 
димость шага 2 вызвана тем, что увеличивающая цепь может проходить по 
дуге несколько раз на разных барьерных уровнях. Благодаря использова- 
нию очереди каждая увеличивающая цепь будет иметь минимальную дли- 
ну. 

Пример 2. Рассмотрим граф с источником 5 и стоком Ё (рис.7), 
пропускные способности всех дуг равны 1. Обозначим поток через дугу 
вектором, первая компонента которого — поток с нулевым барьерным уров- 
нем, вторая — с первым. 


С 
©) 


Рис.7. Начальное состояние графа 


Дуга (6,2) -— барьерная, с барьерным уровнем 1, дуга (6,<) - 
увеличивающая, остальные дуги нейтральные. Применим наш алгоритм. 

После первого шага вершины окажутся помеченными так, как изоб- 
ражено на рис.8. 
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С 
[© 
Н 4 
ор >о >07 
Н Я н Ь в 


Рис. 8. Состояние графа после первого шага 


Таким образом, найдена увеличивающая цепь из источника в сток, 
проходящая по вершинам 5,<,6,с.а,Б,#. При этом по дуге (а,Б) 
цепь проходит дважды: сначала на первом барьерном уровне, затем на 
втором. Поэтому на втором шаге алгоритма величина, на которую по этой 
цепи можно увеличить поток, будет определена как 0,5. После увеличе- 
ния потока и повторения первого шага алгоритма граф будет выглядеть 
следующим образом (рис. 9). 

С 


О 


Н Ян Ь Б 


Рис. 9. Состояние графа, при котором алгоритм завершает работу 


Суммарный поток по дуге (а,Б) равен 1, и эта дуга насыщена. 
Вершина С не может быть помечена, так как вершина ( помечена на ну- 
левом уровне, а поток по дуге (с,<) идет по первому уровню. Пометка 
вершины С и построение через нее увеличивающей цепи означали бы на- 
рушение условий (2). Сток не может быть помечен, и поэтому алгоритм за- 
канчивает свою работу. На основе примера 2 поток, получаемый предло- 
женным алгоритмом, совпадает с потоком, получаемым симплекс-методом 
для системы (3). 
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Го ($. а) = О 
Го Са. Ъ) = О 
Г, Са. Б) = оО 
Го ©, с) = О 
Г, ©, с) = О 
Го (с.а) = О 
Г, Сс.а) = О 
Го Б. 9 = О 
Е $. 90 =оО 
Го ($. а) =1 
Го Са, Б)+ ЕЁ, (а, Ъ) =1 
Го (©, с) + Е, (Б, ©) = 1 
Тоба Ес ау =1 
Го ($. ОЕ, ($, © =1 
Го ($. а) + Ро (с.а) = Ко (а, Ъ) 
Е, (с, а) = Е, (а, Ъ) 
Го (а, Б) = Ро <Б, с) 
Е, (а, Ъ) = ЕЁ, (6, с) + Е, (6, О 
Го (с ьа) = О 
Ро (Б, с) + Е, ($, с) = Е, (с, а) 
Го (5. а) —> пзах 
(3) 
Пример 3. Несмотря на то, что предлагаемый нами алгоритм мо- 
жет находить максимальный поток в сетях с возможной блокировкой, как в 


примере 1, существуют сети, для которых он не находит максимальное ре- 
шение. В этом случае рассмотрим граф, изображенный на рис. 10. 
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С 
0 се ПЬ 
1/3 1/3 
1/ № 


$О 





1/3 О 11/3 


1/2 Е 1/2 


Рис.10. Граф с блокирующим потоком 


Пропускные способности всех дуг равны 1. Дуги (6,<), (е, Л), 
(=2,й) и (.Р?) - увеличивающие, (6,г) и (е,г) - барьерные, с ба- 
рьерным уровнем 2, остальные дуги нейтральные. Поток по не помеченным 


5 2 
на рис.10 дугам считаем равным 0, общий поток на рис.10 равен г 


Установлено, что найденный симплекс-методом максимальный по- 
ток, равен 1. Поток образуется, если все дуги пути, проходящего по верши- 
нам 5,4,е, =,й,А,/, р, а,Ь,, насыщены, а остальные дуги не нагру- 
жены. Так как на рис. 10 дуги (4,е) и (а,Б) насыщены, а (с, а), 
(1,4) и (р,а) не нагружены на нулевом барьерном уровне, то не су- 
ществует не только увеличивающего пути, но и увеличивающей цепи из 
источника в сток. Следовательно, предложенный алгоритм не может найти 


максимальный поток на этом графе. 
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Выводы. Предложенный нами алгоритм основан на идее Форда-Фалкерсо- 
на нахождения увеличивающих цепей. Ранее нами было показано, что для 
поставленной задачи с барьерными ограничениями существуют блокирую- 
щие потоки — потоки, которые не максимальны, но для них не существует 
увеличивающей цепи. Из примера 3 видно, что алгоритм, основанный на 
поиске увеличивающей цепи, не может найти максимальный поток в сети 
на рис. 10. Отличие данной блокировки от блокировки, рассмотренной в 
[2], состоит в том, что блокировка в [2] представлена как отсутствие уве- 
личивающего пути на графе без дополнительных ограничений на пути ча- 
стиц. Тем не менее, решать эту задачу возможно, применяя симплекс-ме- 
тод к системе уравнений и ограничений (2). 

Автор выражает благодарность профессору Я.М.Ерусалимскому за 
постановку задачи. 
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